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ABSTRACT 

This  paper  presents  a  methodology  to  obtain  polarization  profiles  using  a  COTS  (commercial- 
off-the-shelf)  digital  camera.  It  contains  calibration  techniques  to  insure  robust  and  accurate 
determination  of  the  Stokes  parameters  and  the  system  polarization  state.  It  calculates  the  Stokes 
parameters  for  each  scene  pixel  along  with  the  degree  of  polarization,  azimuth  and  ellipticity 
angles  for  each  of  the  RGB  channels.  A  detailed  calculation  shows  that  the  quasi- 
monochromatic  approximation  for  the  wavelength  dependence  of  the  quarter  wave  plate  can  be 
extended  to  include  the  broad  band  RGB  channels.  A  simple  and  inexpensive  camera  calibration 
converts  RGB  values  to  optical  density  and  then  to  relative  intensity  values.  An  analysis  is 
presented  for  the  wavelength  dependent  response  of  the  tuned  zero-order  quartz  quarter-wave 
retarder,  the  broad  band  RGB  channels  and  the  fourth  Stokes  parameter.  The  noise 
characteristics  of  the  digital  camera  are  analyzed  along  with  a  detailed  calibration  of  the  imagery 
for  target  characterization  applications.  An  extensive  validation  process  using  white  light 
compares  known  polarization  states  with  the  empirical  results  from  the  COTS  camera.  A  RGB 
pseudo-coloring  scheme,  related  to  the  PoincarS  sphere,  is  introduced  for  the  visualization  of 
polarization  parameters  of  daylight  scenes. 
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INTRODUCTION 


Our  entry  into  the  polarization  scene  began  over  a  decade  ago.  At  that  time,  polarization 
measurements  were  acquired  using  spot  measurements  from  a  luminance  meter  and  a  typical 
SLR  35  mm  camera1'2.  During  the  previous  four  years,  we  have  studied  the  problems  associated 
with  using  a  COTS  (commercial-off-the-shelf)  digital  camera  to  determine  polarization 
parameters  in  a  daylight  scene3'5.  Although  a  number  of  registration  problems  have  yet  to  be 
completely  resolved,  we  will  show  that  it  is  possible  to  make  polarization  measurements  using  a 
typically  “quasi-monochromatic”  methodology  with  a  COTS  digital  camera.  This  paper 
describes  a  methodology  for  measuring  the  four  Stokes  parameters  in  each  of  the  approximately 
100  nm  RGB  channels  using  a  COTS  digital  camera.  A  calibration  scheme  converts  RGB  pixel 
values  to  relative  intensities  and  the  spatial,  thermal  and  electrical  noise  characteristics  of  the 
CCD  detector  are  examined.  An  equation  compensates  for  the  phase  difference  between  the 
incident  wave  and  the  tuned  frequency  of  a  quartz  zero-order  quarter-wave  retarder.  An 
extensive  validation  process  using  white  light  compares  known  polarization  states  with  the 
empirical  results  from  the  COTS  camera.  An  RGB  pseudo-coloring  scheme,  related  to  the 
Pomcare  sphere,  is  introduced  for  the  visualization  of  polarization  parameters  of  daylight  scenes. 


STOKES  PARAMETER  FORMALSIM 

The  human  visual  system  cannot  easily  detect  polarized  light.  Some  observers  can  determine, 
through  sighting  Haidinger’s  brush6,  that  skylight  is  polarized.  However,  the  human  visual 
system  is  incapable  of  determining  the  complete  polarization  state  of  a  light  beam.  Three 
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independent  parameters  must  be  determined  for  each  scene  pixel  to  completely  characterize  a 
polarization  image7 .  The  Stokes’  parameter  representation  is  one  methodology  for  measuring 
the  polarization  state  of  an  image* .  This  method  involves  four  intensity  measurements  for  each 
scene  pixel  using  four  distinct  optical  configurations.  Each  measurement  corresponds  to  the 
intensity  of  the  beam  after  it  passes  through  each  of  four  different  filter  system  arrangements. 
The  four  Stokes  parameters,  sometimes  called  So,  Si,  S2  and  S3,  are  derived  from  these  four 
measured  intensities  and  form  a  four-element  column  vector  in  four-dimensional  mathematical 
space.  The  Stokes  parameters  are  applicable  to  a  beam  of  light  that  is  completely  polarized, 
partially  polarized  or  unpolarized;  the  beam  may  be  monochromatic  or  polychromatic.  The 
Stokes  parameters  for  completely  polarized  light  propagating  along  the  +z-axis  are 

So  =  E0x2  +  Eoy2  Si  =  E0x2-Eoy2  S2  =  2  Eox  Eoy  cos  8  S3  =  2  Eox  E0Y  sin  5  (1) 

Normalized  parameters  are  obtained  by  dividing  S0,  Si,  S2  and  S3  by  S0.  For  partially  polarized 
light  the  degree  of  polarization  P  is  given  by9 


p  1  Polarized  _  V5.'  +  S2*  +  S^_ 

I  road  “  *0 

where  Ip0iarizcd  is  the  intensity  of  the  polarized  component  and  Ix0tai  is  the  total  intensity.  When 
the  Stokes  parameters  are  normalized,  the  degree  of  polarization  P  becomes  the  radius  of  a 
sphere  with  center  coincident  to  the  typical  Poincare  sphere  (see  Figure  1)  of  unit  radius  (P  1). 
The  inside  of  the  Poincare  sphere  can  be  used  to  represent  the  polarized  portion  of  partially 
polarized  light  (P  <  1).  The  polarization  forms  on  the  surface  of  the  inside  spheres  (P  <  1)  are 
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exactly  the  same  as  the  polarization  forms  on  the  surface  of  the  sphere  of  unit  radius.  From  the 
geometry  of  Figure  1  and  Eq.  (2)  it  can  be  shown  that  for  normalized  Stokes  parameters  and  for 

any  degree  of  polarization 


Si  =  P  cos  2 vp  cos  2x  S2  =  P  sin  2 vj;  cos  2*  S3  =  P  sin  2x 


(3) 


In  fact,  Eq.  (3)  represents  the  spherical  coordinates  of  any  point  on  the  surface  or  inside  of  the 
Poincare  sphere  where  (x,  y,  z)  =  (Si,  S2,  S3). 

The  polarization  form  associated  with  elliptical  polarization,  as  shown  in  Fig.  2,  is  given  by  the 
following  expression 


M 

JE>) 

KEoy  , 

V^o*  > 

-2l 


EyEx 
V  E0yE0x  j 


cos  8  -  sin 2  8 


(4) 


The  polarization  azimuth  angle,  vp,  and  the  ellipticity  angle,  %,  are  also  defined  in  Fig.  2.  Using 


Eqs.  (2)  and  (3),  vp  and  x  can  be  determined  from 


sin  2X  =  -,  =  S\  —  ^  tan2?p  =  ^- 

-x/Sj2  +S22  +S32  (5) 

The  parameter  x  varies  from  +45°  to  -45°;  it  is  positive  for  right-handed  polarization  forms  and 

negative  for  left-handed  polarization  forms.  The  parameter  vp  varies  from  0°  to  180°;  it  is  0°  for 

horizontal  polarization  forms  and  90°  for  vertical  polarization  forms. 
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CAMERA  CALIBRATION 


COTS  digital  cameras  are  effective  “black  boxes”  with  little  or  no  support  from  the  manufactures 
with  regard  to  their  calibrated  response  functions  for  the  various  RGB  channels.  They  do  not 
provide  schematics  for  electronic  or  digital  circuits  to  implement  noise  and  color  correction 
algorithms.  The  COTS  camera  manufacturer  has  different  goals  and  objectives  for  camera 
performance  than  does  the  scientist  or  engineer  who  calculates  the  Stokes  parameters  from 
calibrated  digital  imagery.  The  calibration  process  is  particularly  important  for  extracting  a 
linear  intensity  relationship  from  RGB  grey  level  image  data.  An  accurate  calibration  in  terms  of 
the  intensity  values  insures  that  linear,  spatial  averaging  over  multiple  pixels  improves  signal  to 
noise  characteristics  and  reduces  random  noise  in  the  polarization  imagery. 

Dark  Frames  and  Flat-Field  Frames 

COTS  digital  cameras  typically  use  rectangular  arrays  of  photoreceptors  in  a  sensor 
configuration  called  charge  coupled  devices  (CCDs)10.  Each  CCD  image  consists  of  individual 
detector  elements  or  pixels,  which  accumulate  photoelectric  charges  that  integrate  over  time. 

The  detection  and  amplification  stages  introduce  noise  into  the  imagery  which  consists  of 
thermal  or  “dark  currents”  and  electronic  noise.  An  analysis  of  imagery  taken  with  a  lens  cap 
over  the  camera  lens,  a  dark  image  frame,  is  used  to  quantify  the  effects  of  these  types  of  noise 
phenomena.  The  characteristics  of  a  dark  image  frame  are  different  for  each  kind  of  digital 
camera  and  for  each  of  its  RGB  channels.  Ideally,  the  non-zero  values  obtained  from  a  dark 
frame  need  to  be  subtracted  out  of  the  non-dark  image  frames.  An  analysis  of  the  dark  image 
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frame  of  an  Epson  850Z  digital  camera  showed  absolutely  no  non-zero  values  for  any  of  its 
channels.  For  this  camera,  thermal  noise  is  subtracted  out  prior  to  the  determination  of  the  RGB 
values.  Typically,  non-zero  values  obtained  from  a  dark  frame  change  with  exposure  time  and 
zoom  setting.  A  comprehensive  study  of  this  phenomenon  is  essential  for  accurate  scientific 
work. 

Flat  frames  are  CCD  images  created  by  exposing  the  CCD  to  a  uniform  field  of  light.  They  can 
be  used  to  calibrate  variations  in  photosite  sensitivity  for  the  entire  CCD  chip.  These  images  are 
a  photosite  map  of  the  CCD’s  sensitivity  to  light.  If  large  variations  in  sensitivity  occur 
consistently  at  a  particular  photosite,  they  can  be  corrected  with  a  computer  program. 

Spectral  Response  of  the  RGB  Channels 

A  study  was  conducted  to  determine  the  spectral  response  of  the  RGB  channels  of  the  relatively 
inexpensive  Epson  850Z  digital  camera.  Twelve  different  narrow  band  (10  nm)  color  filters 
(diameter  =  2  inches)  ranging  from  440  nm  to  660  nm  were  positioned  sequentially  in  front  of 
the  camera.  A  collimated  beam  of  white  light  was  transmitted  through  the  filters  and  into  the 
camera  producing  a  uniformly  colored  circular  image.  Since  adjacent  pixels  have  slightly 
different  sensitivities,  average  RGB  values  were  calculated  using  a  10  X  10  pixel  area  from  each 
digital  image.  This  procedure  was  followed  for  each  of  the  color  filters  used  in  the  study.  Figure 
3  gives  an  example  of  the  average  RGB  values  obtained  from  each  of  the  three  channels  when  in 
the  fixed  white  balance  (W/B)  mode  of  the  Epson  850Z. 


6 


As  shown  in  Figure  3,  the  red  channel  response  has  a  sharp  cut-off  at  approximately  600  nm. 

The  red  channel  continues  to  respond  to  wavelengths  up  to  approximately  540  nm;  however,  a 
small  response  is  observed  for  wavelengths  less  than  480  nm  for  this  channel.  Figure  3  shows 
that  the  green  channel  responds  to  all  wavelengths  between  440  nm  and  660  nm  with  a  peak 
response  at  540  nm.  The  green  channel  is  most  sensitive  to  wavelengths  between  520  nm  and 
560  nm.  The  blue  channel  response  has  a  sharp  cut-off  at  approximately  480  nm.  The  blue 
channel  continues  to  respond  to  wavelengths  up  to  approximately  520  nm;  however,  a  small 
response  is  observed  for  wavelengths  greater  than  600  nm  for  this  channel. 

RGB  to  Optical  Density  to  Relative  Intensity 

The  digital  camera  converts  each  8-bit  RGB  channel  value  into  an  integer  between  0  and  255. 
These  RGB  pixel  values  are  not  intensity  values  that  can  be  used  to  immediately  calculate  the 
Stokes  parameters.  To  determine  the  mathematical  nature  of  the  conversion  of  RGB  pixel  values 
to  light  intensity,  a  Macbeth  Color  Checker1 1  was  photographed  with  different  neutral  density 
filters  in  front  of  the  lens  of  the  digital  camera.  Neutral  density  filters  have  known  optical 
densities,  which  reduce  the  transmission  through  them  according  to 

It  =  Io  10°  (5) 

where  D  =  optical  density,  I0  =  incident  intensity  and  IT  =  transmitted  intensity.  The  achromatic 
gray  level  portions  of  the  Macbeth  Color  Checker  (colors  #19  -  #24)  have  a  broadband  spectral 
response  with  reflectance  values  corresponding  to  known  optical  densities.  Figure  4  shows  the 
average  green  channel  pixel  values  for  different  product  combinations  of  Macbeth  reflectance 
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and  neutral  density  filter  transmission.  A  MatLab12  algorithm  calculated  the  average  pixel 
response  over  a  10  X  10  pixel  area.  The  best  polynomial  curve  fit  through  the  points  was 
determined  to  be  a  hyperbola  of  the  form 

Ax2  +  Bxy  +  Cy2  +  Dx  +  Ey  +  F  =  0  (6) 

where  x  and  y  correspond  to  optical  densities  and  RGB  values,  respectively.  The  coefficients  in 
Eq.  (19)  are  different  for  each  color  channel  of  the  digital  camera.  Using  the  fixed  white  balance 
option  of  the  Epson  850Z,  the  coefficients  obtained  for  the  green  channel  are: 

A  =  C  =  0.19926,  B  =  35.67852,  D  =  127.49000,  E  =  -55.65087  and  F  =  -520.1 1 10  (7) 

MEASURING  THE  STOKES  PARAMETERS 


The  Stokes  parameters  can  be  measured  using  a  linear  polarizer  and  a  retarder.  If  light  is  first 
transmitted  through  a  linear  retarder  with  its  fast  axis  oriented  at  an  angle  O  with  respect  to  the 
x-axis  and  then  transmitted  through  a  linear  polarizer  with  its  transmission  axis  oriented  at  an 
angle  9  with  respect  to  the  x-axis,  the  intensity  of  the  transmitted  light  can  be  expressed  as 

I(Q,0,e)  =  0.5  {So  +  [cos  29  -  sin  2Q  sin  2(Q  -  9)  (1  -  cos  s)]  Si 

+  [sin  2Q  cos  2(Q  -  0)  -  cos  20  sin  2(0  -  9)  cos  e  ]  S2  -  sin  2(0.  -  0)  sin  8  S3}  (8) 
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where  e  is  the  phase  difference  produced  by  the  retarder  and  is  a  function  of  the  retarder 
thickness.  I(Q,9,s)  denotes  an  intensity  measurement  corresponding  to  a  particular  set  of  values 
for  Q,  0  and  e.  Equation  9  compensates  for  the  difference  in  phase  shifts  between  the  incident 
light  and  a  retarder  phase  shift  of  90°  where  X  is  the  wavelength  of  the  incident  light  and  Xj  the 
tuned  wavelength  of  the  retarder,  both  in  nanometers1 1 . 


Aj.  -50.876^ 
2  ^  /l  -50.876 , 


(9) 


Figure  5  shows  the  phase  difference  for  a  quartz  retarder  as  a  function  of  the  input  wavelength 
for  specific  tuned  wavelengths.  The  slope  ds/dX  in  Figure  5  has  its  smallest  values  for  a  quartz 
retarder  tuned  to  the  blue  end  of  the  spectrum  and  for  input  wavelengths  in  the  red  end  of  the 
visible  spectrum  (the  red  channel  of  a  COTS  digital  camera). 


The  Stokes  parameters  can  be  determined  experimentally  from  the  following  four  measurements, 
as  shown  in  Figure  6: 


Ii=I(0,0,0)  h  ~  1(0,  90, 0)  I3  =  1(0, 45,0)  I4  =  I(0,45,e)  (10) 

In  this  arrangement,  a  retarder  is  inserted  into  the  optical  path  for  only  the  fourth  measurement. 
The  four  measurements  described  in  Eq.  (10)  are  not  the  only  four  measurements  that  will  lead  to 
a  determination  of  the  Stokes  parameter.  In  a  previous  work13,  we  describe  a  simple  two- 
element  variable  filter  system  where  neither  the  polarizer  nor  the  retarder  needs  to  be  removed 
for  any  measurement. 
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Figure  7  shows  an  experimental  setup  to  acquire  this  four-image  sequence.  It  shows  a  digital 
camera,  a  linear  polarizer,  a  hinged  retarder  and  a  hinged  filter  mounted  on  a  tripod.  The  retarder 
and  filter  are  hinged  so  that  they  can  be  moved  quickly,  precisely  and  reproducibly  into  and  out 
of  the  optical  path  perpendicular  to  the  optical  axis.  Substituting  Eq.  (10)  into  Eq.  (8)  confers 
that 

50  =  1(0,  0,  0)  + 1(0,  90,  0)  =  Ii  +  I2 

51  =  1(0, 0, 0)- 1(0, 90, 0)  =  I,  - 12  ; 

52  =  2  1(0,  45,  0)  -  S0  =  2  I3  -  S0 

c  _  2/4  -S0  -S2coss  (11) 

sm  £ 

The  following  sample  calculations  determine  the  degree  of  polarization,  P,  the  polarization 
azimuth  angle,  vp,  and  the  polarization  ellipticity  angle,  x>  for  a  single  pixel  in  the  green  channel 
of  the  Epson  850Z  digital  camera.  Assume  the  following  data:  (1)  the  tuned  wavelength  for  the 
retarder  is  436  nm,  (2)  the  input  wavelength  is  550  nm  and  (3)  the  RGB  pixel  values 
corresponding  to  Ij,  I2, 13  and  I4  are  125,  123,  137  and  127  respectively.  Using  Eq.  (9)  gives  £  = 
1.212  radians.  Using  Eqs.  (6)  and  (7)  to  convert  RGB  green  channel  values  to  optical  density 
yields  0.951,  0.963,  0.878  and  0.929.  The  corresponding  relative  intensities,  using  Eq.  (5),  are: 
Ii  =  0.1 12, 12  =  0.109, 13  =  0.132  and  I4  =  0.1 18.  Calculations  of  the  Stokes  parameters  from  Eq. 
(11)  yields  (So,  Si,  S2,  S3)  =  (0.221,  0.003,  0.044,  -0.001).  The  normalized  Stokes  parameters  are 
(S0*,  Si*,  S2*,  S3*)  =  (1.000,  0.014,  0.200,  -0.003).  From  Eqs.  (2)  and  (3):  (P,  \\ /,  x)  =  (0.200, 43 
deg.,  -0.4  deg.). 
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BROADBAND  RELATED  ERRORS 


The  phase  difference  of  a  retarder  is  wavelength  dependent,  as  Eq.  (9)  illustrates.  Hence,  the 
fourth  Stokes  parameter,  expressed  in  Eq.  (11),  is  also  wavelength  dependent.  Since  Eq.  (9)  can 
only  accommodate  one  input  wavelength  and  the  bandwidth  of  an  RGB  channel  is  approximately 
100  nm,  errors  are  introduced  into  polarization  parameters  dependent  on  S3;  namely,  P  and  % 
(Eqs.  2  and  5). 


Although  s  and  A  are  not  related  linearly  throughout  the  visible  spectrum,  they  are  nearly  so 
within  the  approximate  100  nm  bandwidth  of  an  RGB  channel  (see  Figure  5).  Therefore,  to  a 
good  approximation,  the  wavelength  corresponding  to  the  average  value  of  8  for  a  given  RGB 
channel  is  the  one  corresponding  to  the  midpoint  of  the  channel.  Alternately,  the  average  value 
of  s,  using  a  quartz  quarter-wave  retarder,  can  be  obtained  from 


1  fatt  / l,.  -  50.876 
~  VA  J L  2  A  -50.876  dA 


(12) 


where  Aa  and  Ab  are  the  beginning  and  ending  wavelengths  of  an  RGB  channel.  Integration  of 
Eq.  (12)  gives 


£avg 


v  2  > 


/L  -  50.876  r,  ..  .  ,  ,  1 

- - Lin  (4  -  50.876)  -  In  (Aa  -  50. 876) J 


(13) 


Ab  -  A 


The  input  wavelength  that  corresponds  to  eavg,  Ac,  is  obtained  by  substituting  8avg  for  8  in  Eq.  (9) 
and  solving  for  A,.  Interestingly,  the  value  obtained  for  Ac  is  independent  of  the  tuned 
wavelength  At-  For  Aa  =  400  nm  and  Ab  =  500  nm  (the  blue  channel)  Ac  =  447.9  nm;  for  Aa  = 
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500  nm  and  /.r  =  600  nm  (the  green  channel)  Xc  —  548.3  nm;  for  XA  =  600  nm  and  a.b  =  700  nm 
(the  red  channel)  Xc  =  648.6  nm.  In  other  words,  the  wavelength  corresponding  to  savg  for  a 
given  RGB  channel  is  the  midpoint  wavelength. 

The  fourth  Stokes  parameter,  as  expressed  in  Eq.  (11),  is  the  only  Stokes  parameter  that  is 
wavelength  dependent.  Therefore,  it  is  responsible  for  any  errors  resulting  from  extending  quasi- 
monochromatic  theory  to  include  wider  bandwidths.  The  source  of  bandwidth  errors  originates 
with  I4  in  Eq.  (11).  Solving  Eq.  (1 1)  for  I4  gives 

I4  =  0.5  (So  +cos  s  S2  +  sin  s  S3)  (14) 


Suppose  light  consisting  of  a  flat  spectral  response  for  each  RGB  channel  and  having  the  same 
Stokes  parameters  for  each  wavelength  is  transmitted  though  a  quartz  quarter-wave  retarder.  The 
average  value  of  Lt  is  obtained  from 


I4avg  =  °'5  ~  (S0  +  cos  c  +  sin  e  S3 )  de 


(15) 


where  sA  and  Sb  are  the  phase  differences  corresponding  to  the  beginning  and  ending 
wavelengths  of  an  RGB  channel.  Integration  of  Eq.  (15)  gives 


0.5 


*4avg 


-  £a 


[(sB  -  £  A  )  S0  +  (sin  SB-sin£A)S2  -  (cos  £B  -  cos  SA  )  S3  ] 


(16) 


Table  1  provides  three  examples  of  the  errors  introduced  in  the  fourth  Stokes  parameter  by 
transmitting  a  broadband  polarized  light  beam  through  the  fourth  filter  of  Figure  6.  The  best 
case  and  worst  case  scenario  are  given  for  each  example.  The  first  two  rows  refer  to  Example  1, 
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rows  three  and  four  refer  to  Example  2  and  rows  five  and  six  refer  to  Example  3.  Example  1 
relates  to  right  circular  polarization  with  (S0,  Si,  S2,  S3)  =  (1, 0,0,1)  and  (P,  x)  =  0,  45  deg.). 
Example  2  relates  to  elliptical  polarization  with  (So,  Si,  S2,  S3)  =  (1,-0.250,  0.433,  0.866)  and  (P, 
X)  =  (1,  30  deg.).  Example  (3)  relates  to  linear  polarization  with  (So,  Si,  S2,  S3)  =  (1,-0.500, 
0.866,  0.000)  and  (P,  x)  =  (1,  0  deg.).  The  best  case  scenario  uses  a  %  wave  retarder  tuned  to 
400  nm  and  corrected  to  the  center  frequency  of  the  red  channel.  The  worst  case  scenario  uses  a 
similar  tuned  wavelength  of  700  nm  corrected  to  the  center  frequency  of  the  blue  channel.  In 
Example  1,  the  best  case  scenario  produces  an  error  in  the  value  of  P  that  is  less  than  0.5%  and 
an  error  in  the  value  of  x  that  is  less  than  0.1  degrees;  the  worst  case  scenario  produces  a  7.8% 
error  in  the  P-value  and  an  error  in  the  value  of  x  that  is  less  than  0.1  degrees.  In  Example  2,  the 
best  case  scenario  produces  an  error  in  the  value  of  P  that  is  less  than  0.01%  and  an  error  in  the 
value  of  x  that  is  less  than  0.1  degrees;  the  worst  case  scenario  produces  a  6.3%  error  in  the  P- 
value  and  a  l.l  degree  error  in  the  value  of  X-  In  Example  3,  the  best  case  scenario  produces  an 
error  in  the  value  of  P  that  is  less  than  0.01%  and  a  0.2  degree  error  in  the  value  of  x;  the  worst 
case  scenario  produces  an  error  in  the  value  of  P  that  is  less  than  0.01%  and  a  0.3  degree  error  in 
the  value  of  %. 

EXPERIMENTAL  VALIDATION 


The  approximate  bandwidth  of  white  light  is  between  400  nm  and  700  nm.  Several  COTS 
digital  cameras  were  examined  and  found  to  have  RGB  channels  with  approximate  half-widths 
of  50  nm.  The  typical  spectral  response  of  these  channels  is  approximately  as  follows:  red 
channel:  600-700  nm,  green  channel:  500-600  nm  and  blue  channel:  400  and  500  nm.  Each 
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wavelength  in  a  white  light  beam  may  have  a  different  set  of  Stokes  parameters  and  therefore  a 
different  polarization  form  for  each  narrow  band  region  of  the  spectrum.  Although  there  may  be 
mechanisms  that  produce  abrupt  changes  or  discontinuities  in  the  Stokes  parameters  at  a  specific 
wavelength,  many  processes  produce  uniform  or  slowly  varying  changes  in  the  polarization 
parameters  within  a  white  light  beam.  For  example,  the  normalized  Stokes  parameters  for 
monochromatic  light  transmitted  through  a  linear  polarizer  and  retarder  (shown  in  Figure  8)  are 
given  by 


50  =  1 

51  =  cos  2((3  -  a)  cos  2p  +  sin  2(P  -  a)  sin  2P  cos  a 

52  =  cos  2(p  -  a)  sin  2p  -  sin  2(p  -  a)  cos  2P  cos  a 

53  =  sin  2(p  -  a)  sin  a  (17) 

where  a  is  the  phase  difference  produced  by  the  retarder.  For  a  quartz  retarder1 1 


7rf\  -50.876' 
2^/1-50.876, 


(18) 


Column  A  of  Figure  9  contains  the  superimposed  polarization  forms  associated  with  the 
minimum,  maximum  and  average  wavelengths  within  the  RGB  channels  for  white  light 
transmitted  through  the  apparatus  in  Figure  8.  These  theoretical  values  are  compared  with  the 
experimental  data  in  Column  B  to  validate  the  results  from  the  Epson  850Z  camera.  The 
imagery  was  recorded  using  the  filter  system  of  Figure  6  and  the  apparatus  in  Figure  7.  The 
quartz  retarder  was  a  zero-order,  quarter-wave  plate  tuned  to  546  nm.  Images  were  acquired  for 
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different  elliptically  polarized  states  generated  with  the  apparatus  of  Figure  8  for  each  of  the 
RGB  channels.  A  MatLab  script  converted  RGB  values  to  intensity  for  each  pixel.  The  average 
Stokes  parameters  were  calculated  over  a  10  X  10  pixel  region  to  improve  spatial  signal  to  noise 
characteristics. 

Column  B  compares  the  theoretical  polarization  forms,  using  the  mean  wavelength  in  a  RGB 
channel,  to  empirical  values  obtained  from  the  Epson  850Z.  The  graphical  comparisons  in 
column  B  show  good  agreement  between  the  elliptical  polarization  states  for  ‘quasi- 
monochromatic”  and  broadband  white  light. 

VISUALIZING  POLARIZATION  PARAMETERS  IN  A  DAYLIGHT  SCENE 

The  Stokes  parameters  can  easily  be  encoded  in  a  daylight  scene  by  assigning  RGB  values  to  the 
normalized  values  of  Si,  S2  and  S3  at  each  pixel  site  in  the  scene  as  follows: 

R  =  int[127.5  (1  -  Si)],  G  =  int[127.5  (1  -  S2)]  and  B  =  int[127.5  (1  -  S3)]  (19) 

This  pseudo-color  scheme  closely  relates  to  the  Poincare  sphere  representation  of  polarized  light 
where  each  Stokes  vector  maps  into  unique  RGB  values  for  0*P*  1 .  The  surface  of  the  sphere 
consists  of  totally  polarized  light  with  P  =  1  (Figure  10a).  Unpolarized  light  (Si  =  S2  =  S3  =  0) 
corresponds  to  middle  gray  (R  =  G  =  B  =  127)  at  the  center  of  the  sphere  while  the  inside  of  the 
sphere  relates  to  partially  polarized  light  (0<P<1)  (Figure  10b).  In  this  color  scheme, 
unpolarized  or  weakly  polarized  light  is  middle  gray  or  highly  unsaturated  in  the  primary  colors. 
The  azimuth  and  ellipticity  polarization  angles  are  essential  parameters  in  obtaining  a  complete 
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polarization  profile.  One  method  of  displaying  these  calculated  parameters  is  to  assign  a 
different  color  to  each  specific  angle.  The  equator  of  the  Poincare  sphere  plays  a  very  special 
role  in  our  pseudo-coloring  scheme  (Figure  10c).  It  corresponds  to  linearly  polarized  light  (S3  = 
0)  and  is  also  used  to  encode  the  polarization  azimuth  and  ellipticity  angles  in  a  daylight  scene, 
as  given  in  Eq.  (20).  Substituting  x  for  \\)  in  Eq.  (20)  produces  a  color-mapping  scheme  for  the 
X-images  as  given  in  Eq.  (21). 

R  =  int[  127.5  (1  -  cos  2 v|/)],  G  =  int[127.5  (1  -  sin  2\j /)]  and  B  =  127  (20) 

R=  int[127.5  (1  -  cos  2x)],  G  =  int[127.5  (1  -  sin  2x)]  and  B  =  127  (21) 

Since  black  (0,0,0),  white  (255,255,255)  and  yellow  (255,255,0)  are  excluded  from  the  encoding 
scheme  of  Eqs.  (20-21),  they  are  used  for  special  conditions.  Black  is  used  when  S0  =  0,  white 
for  indeterminate  values  such  as  dividing  by  zero  and  for  P  =  0;  yellow  is  used  in  the  ^-images 
whenever  x  =  ±  45°  (circular  polarization).  The  normalized  Stokes  parameters  for  circular 
polarized  light  are  (S0,  Si,  S2,  S3)  =  (1,0,0,±1).  Since  the  y-value  is  related  to  the  ratio  S2/S1  (Eq. 
5),  \\i  is  indeterminate  for  circular  polarization.  Polarization  states  that  are  nearly  circular  are  ill- 
conditioned  with  respect  to  y.  Referring  to  Eq.  (5),  if  sin  2x  is  greater  than  one  or  less  than  -1 
(due  to  measurement  errors),  then  x  is  undefined.  The  undefined  x-values  are  circumvented  by 
assuming  x  =  45°  when  sin  2x  >  1  and  x  =  -45°  when  sin  2 x  <  -1. 

The  degree  of  polarization,  P,  varies  between  0  and  1 .  The  simplest  method  of  encoding  this 
parameter  in  a  daylight  scene  is  to  use  the  equation  [pixel  value]  =  255  P.  The  black  areas  (pixel 
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value  =  0)  in  the  resulting  monochrome  image  will  correspond  to  no  polarization  and  the  white 
areas  (pixel  value  =  255)  will  correspond  to  light  that  is  100  percent  polarized. 

Figure  1 1  gives  an  example  of  the  use  of  the  visualization  process  described  above.  The 
polarization  parameters  P,  \p  and  x  are  displayed  for  a  vehicle  in  daylight  at  12:00  on  April  26, 
2000  in  Warren,  MI.  The  altitude  of  the  sun  was  55  degrees  and  the  azimuth  (south  =  0)  of  the 
sun  was  319  degrees.  The  interpretation  of  these  polarization  parameters  as  a  function  of  the 
relative  position  of  the  sun  have  been  analyzed  and  compared  to  laboratory  data.  These  findings 
will  be  presented  in  a  future  paper. 

CONCLUSIONS: 


This  paper  develops  a  systematic  methodology  to  measure  the  four  Stokes  parameters  for  each 
pixel  in  an  RGB  image  where  the  scene  is  illuminated  by  broad  band  radiation  in  the  visual  part 
of  the  electromagnetic  spectrum.  A  calibration  procedure  maps  the  COTS  digital  camera  RGB 
eight  bit  grey  levels  into  their  corresponding  relative  intensity  values.  The  calibration  process 
produces  linear  intensity  values  that  allow  for  linear  averaging  over  a  selected  number  of  pixel 
regions  of  the  image.  This  process  reduces  random  errors  in  the  polarization  images.  This 
capability  is  very  important  because  of  the  inherent  ill-conditioning  in  many  polarization  states. 
The  fourth  Stokes  parameter  measurement  requires  a  retarder  element,  which  is  wavelength 
dependent  and  normally  limited  to  narrow  band  incident  light.  Detailed  calculations  show  that 
the  quasi-monochromatic  approximation  can  be  extended  to  include  the  spectral  response  of  the 
broad  band  RGB  channels.  We  provide  equations  to  compensate  for  the  phase  differences 
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between  the  tuned  frequency  of  the  retarder  and  an  arbitrary  input  wavelength.  Excellent 
agreement  is  achieved  between  the  theoretical  and  experimental  data  for  known  elliptical 
polarization  states. 

The  Stokes  parameters  are  encoded  into  a  polarization  image  by  assigning  RBG  color 
coordinates  to  their  normalized  values.  This  color  encoding  scheme  closely  relates  to  the 
Poincare  sphere  representation  of  polarized  light  where  each  Stokes  vector  maps  into  unique 
RGB  values  for  0  <  P  <  1.  We  then  use  this  visualization  technique  to  represent  the  degree  of 
polarization,  azimuth  and  ellipticity  angles  associated  with  natural  daylight  scene  imagery. 

Future  work  will  examine  the  following  issues  arising  from  this  research.  Interpretation  of  the 
physical  phenomena  associated  with  polarization  imagery  including  daylight  illumination 
conditions  and  material  characteristics  within  the  scene.  We  will  analyze  the  effects  of  spatial 
and  temporal  registration  errors  upon  the  fidelity  of  P,  \\i  and  x  polarization  images. 
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Table  Caption 


Table  1.  Three  examples  of  the  errors  introduced  in  the  measurement  of  the  fourth  Stokes 
parameter  by  transmitting  a  broadband  polarized  light  beam  through  the  fourth  filter  of  Figure  6. 
The  first  two  rows  refer  to  Example  1,  rows  three  and  four  refer  to  Example  2  and  rows  five  and 
six  refer  to  Example  3.  The  best  case  and  worst  case  scenario  are  given  for  each  example.  The 
best  case  scenario  uses  a  tuned  wavelength  of  400  nm  and  the  red  channel.  The  worst  case 
scenario  uses  a  tuned  wavelength  of  700  nm  and  the  blue  channel. 

Figure  Captions 

Figure  1.  The  Poincare  sphere. 

Figure  2.  The  polarization  ellipse. 

Figure  3.  The  channel  spectral  response  of  an  Epson  850Z  digital  camera. 

Figure  4.  Conversion  from  RGB  to  optical  density  for  the  green  channel  of  an  Epson  850  Z.  A 
Matlab  algorithm  calculated  the  average  pixel  response  from  a  10  X  10  pixel  area. 

Figure  5.  The  relationship  between  the  phase  difference  of  a  zero-order  quartz  quarter  wave 
retarder  and  the  wavelength  transmitted  through  it  for  four  different  tuned  wavelengths:  650,  550 
and  450  nm. 

Figure  6.  The  four-filter  system  used  to  determine  the  Stokes  parameters.  P  is  the  transmission 
axis  of  a  linear  polarizer  and  R  is  the  fast  axis  of  a  zero  order  quartz  quarter-wave  retarder. 

Figure  7.  An  experimental  camera  setup  to  acquire  the  four-image  sequence  for  obtaining  the 
Stokes  parameters  associated  with  each  pixel  in  the  scene. 
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Figure  8.  Apparatus  for  producing  known  polarization  forms.  Light  of  wavelength  X  is 
transmitted  through  a  linear  polarizer  with  its  transmission  axis  at  an  angle  a  relative  to  the  x- 
axis  and  transmitted  through  a  zero-order  quarter-wave  retarder  with  its  fast  axis  at  an  angle  P  to 
the  x-axis. 

Figure  9.  Column  A:  Polarization  forms  generated  with  the  apparatus  of  Figure  8  and 
corresponding  to  the  minimum,  maximum  and  midpoint  wavelengths  within  each  color  channel. 
Column  B:  A  comparison  between  the  known  polarization  forms,  produced  by  using  a  midpoint 
wavelength  for  each  RGB  channel,  to  those  actually  obtained  from  the  Epson  850Z. 

Figure  10.  (a)  The  surface  of  the  colorized  Poincare  sphere,  (b)  The  colorized  Poincare  sphere 
of  unit  radius.  On  the  surface  of  the  sphere,  P  =  1.  At  the  center,  P  =  0.  Inside  the  sphere  P  <  1. 
(c)  The  Equatorial  plane  of  the  colorized  Poincare  sphere. 

Figure  11.  An  example  of  the  use  of  the  encoding  methodology  to  colorize  the  polarization 
parameters.  The  degree  of  polarization  P,  the  polarization  azimuth  angle  \\i  and  the  ellipticity  y_ 
are  shown  for  12:00.  Viewing  is  toward  the  north. 
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Channel  Spectal  Response 
for  the  Epson  850Z 
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Figure  7 
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Figure  8 
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Figure  9 
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(c)  x  at  12:00 


Figure  10 


Figure  11 


